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ABSTRACT 

Three-dimensional stellar modeling has enabled us to identify a deep-mixing 
mechanism that must operate in all low mass giants. This mixing process 
is not optional, and is driven by a molecular weight inversion created by the 
3 He( 3 He,2p) 4 He reaction. In this paper we characterize the behavior of this mix- 
ing, and study its impact on the envelope abundances. It not only eliminates the 
problem of 3 He overproduction, reconciling stellar and big bang nucleosynthesis 
with observations, but solves the discrepancy between observed and calculated 
CNO isotope ratios in low mass giants, a problem of more than 3 decades' stand- 
ing. This mixing mechanism, which we call '<5/i-mixing', operates rapidly (relative 
to the nuclear timescale of overall evolution, ~ 10 8 yrs) once the hydrogen burn- 
ing shell approaches the material homogenized by the surface convection zone. 
In agreement with observations, Pop I stars between 0.8 and 2.0 M & develop 
12 C/ 13 C ratios of 14.5 ± 1.5, while Pop II stars process the carbon to ratios of 
4.0 ± 0.5. In stars less than 1.25 M & , this mechanism also destroys 90% to 95% 
of the 3 He produced on the main sequence. 

Subject headings: stars: red giants; abundance anomalies 



- 2 - 



1. Introduction 

In a previous paper (Dearborn et al. 2006; hereinafter Paper I) we used a fully 3- 
dimensional code Djehuty, developed in the Lawrence Livermore National Laboratory, to 
investigate the onset of the helium flash in a low- mass red giant. While convective motion due 
to the helium flash was seen to occur in much the same region as predicted by one-dimensional 
(spherically symmetric, hydrostatic) models, we noticed some minor motion, apparently also 
of a turbulent convective character, in an unexpected region: a region above but not far above 
the hydrogen-burning shell, and well below the base of the conventional surface convection 
zone. This is visible in Fig. 14 of Paper I. On subsequent close inspection we found that 
this additional motion was due to a very small molecular-weight inversion which developed 
into a dynamical instability that we identified as a Rayleigh- Taylor instability (Eggleton et 
al. 2006, 2007; hereinafter Papers II, III). This inversion was due to the burning of 3 He, of 
which quite a high concentration is left by the retreating surface convection zone. 

In Papers II and III we showed why such an inversion is expected to arise, once the 
surface convection zone, having reached its deepest extent (the 'First Dredge-Up', or FDU), 
begins to retreat. Furthermore, we suggested that this motion should grow in extent so 
that it reaches upwards from the zone where the //-inversion is maximal (which is also the 
zone where 3 He-burning is maximal) to the base of the normal surface convection zone. It 
should lead to the destruction of most (> 90%) of the 3 He in the surface layers, and it 
should simultaneously allow for some processing of 12 C to 13 C. Thus at the surface of the 
star 3 He should be progressively depleted, and 13 C progressively enhanced, beyond the values 
expected from conventional ID models. 

In this paper we consider further the effect of our additional mixing, which we refer to 
as '<5//-mixing'. We would like to emphasise that this mixing can explain in a natural way 
the observed abundances that have hitherto been attributed to mechanisms like rotation and 
magnetic fields. We would also like to emphasise the value of 3D modeling. It was only by 
using a fully 3D (hydrodynamic) model that this mixing was noticed. It is due to a very 
slight effect, a //-inversion that amounts to less than one part in 10 4 , and yet it is quite 
obvious in a 3D simulation. One can in fact see the inversion in ID models - although we 
are not aware that anyone has actually commented on it - but in ID models mixing only 
occurs if the code-developer tells the code to include it. 

Three dimensional simulations are expensive of computer time. Our philosophy in using 
D j ehuty has been that 3D simulations spanning a short period of time may allow us insight 
into complex hydrodynamic and hydromagnetic phenomena which will then enable us to 
improve the quality of simplistic ID models. We follow this principle here. We introduce 
into our ID code a convective mixing coefficient that depends on the //-gradient, but only 
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if the /i-gradient is in the sense that it is destabilising. This is in addition to the mixing 
coefficient that is due to ordinary convective instability, which depends (crudely speaking) 
on the entropy gradient but also only if this gradient has the sign which is destabilising. 

In Section 2 we discuss the significance of possible mixing for abundance measurements 
on the Giant Branch. In Section 3 we discuss our new 5/x-mixing mechanism in more detail. 
In Section 4 we describe a ID code which incorporates a simple model of our 5/i-mixing. In 
Sections 5 and 6 we discuss the sensitivity of this model to some input parameters that we 
estimate. In Sections 7 - 12 we present our results and conclusions. 



2. The Importance of Mixing on the Giant Branch 

There is a long history of observing the carbon, nitrogen and oxygen (CNO) isotope 
ratios in Red Giants (Lambert & Dearborn 1972, Day et al. 1974, Dearborn et al. 1975, 
1976, Tomkin et al. 1975, 1976, Tomkin & Lambert 1978, Harris & Lambert 1984a,b, Gilroy 
1989) as a probe of stellar interiors and evolution. While on the main sequence, nuclear 
reactions change the abundance distribution in a star's deep envelope. These changes are 
then mixed to the surface in the FDU, when the star becomes a giant. Among the isotopes 
that are substantially enhanced by low mass stars are 3 He, 13 C and 17 0. Other isotopes like 
15 N and 18 are reduced. Distressingly, the observed ratios were nearly always different from 
what was predicted, and sometimes very different (Dearborn et al. 1976, Dearborn 1992). 

A result of giant-branch abundance changes that is of particular consequence is the 
robust result that low-mass stars (below ~ 1.5 M ) are major producers of 3 He. In these 
stars, the PP chain produces a 3 He-rich peak in the envelope of the star. Giant-branch 
convection homogenizes this region, in an episode known as the First Dredge-Up or FDU, 
raising the surface 3 He abundance by a factor of ~ 80 (see Fig. 1). As the helium core 
grows some 3 He is destroyed, but only (in classical ID models) below the surface convection 
zone. The fractional abundance of 3 He in the convective envelope itself is not diminished, 
even though the mass of the convective envelope is diminished. The presence of horizontal 
branches in the HR diagrams of globular clusters can be understood only if substantial 
amounts of envelope mass are lost on the giant branch. This should result in a substantial 
enhancement of the 3 He in the interstellar medium, by low-mass stars of both Pop I and 
Pop II (Dearborn et al. 1986, Steigman et al. 1986, Dearborn et al. 1996). Work by Hata 
et al. (1995) argued that unless > 90% of this 3 He is destroyed prior to ejection, the results 
of stellar nucleosynthesis come into conflict with standard Big Bang Nucleosynthesis. To 
avoid this, there must be some mechanism that operates well before the helium core flash to 
destroy the 3 He. 
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Sweigart & Mengel (1979), Smith & Tout (1992) and Wasserburg et al. (1995) posited 
that a deep rotationally-driven circulation mechanism could solve the problem of low 12 C/ 13 C 
values observed in red giants, and destroy the excess 3 He at the same time. Mechanisms 
have been explored that might cause this deep mixing, including magnetic-field generation 
(Hubbard & Dearborn 1980) or rotational mixing (Charbonnel 1995, Chaname et al. 2005). 
Note in partricular that Palacios et al. (2006) have shown recently that rotation is unable 
to produce the required mixing. The necessity for some deep mixing was also discussed by 
Hogan (1995). 

To destroy the necessary 3 He, any such mechanism must operate efficiently in nearly 
all low-mass stars. While rotation is certainly present in young stars it is inferred to decay 
on the main sequence. The slow rotation speeds seen in white dwarfs, and the rotation 
rate inferred in the Solar core by helioseismology, suggest that angular momentum is not 
conserved in the cores either. As a consequence, dependence on the ability of rotation to 
destroy the excess 3 He in all long-lived stars is unsettling. Nevertheless, rotational mixing 
must surely exist in some stars, though of poorly known efficiency. 

Here, we wish to examine the implications of our new mixing mechanism, described 
in Papers II and III, that operates efficiently in all low-mass stars. This new mechanism 
is not optional: it inevitably arises on the First Giant Branch when the hydrogen-burning 
shell encroaches on the homogenised, formerly convective, zone left behind by the retreating 
convective envelope, and begins to burn 3 He. This 3 He burning occurs just outside the normal 
hydrogen-burning shell, at the base of a radiatively stable region about 1 R & in thickness. 
The burning causes a molecular weight inversion, which drives mixing all the way to the 
usual convection zone. As we will show, it not only routinely destroys > 90% of the 3 He 
produced in the low-mass offenders, but reduces the 12 C/ 13 C ratios from the range of 20 to 
35 down to the observed values between 5 and 15, depending on metallicity. 



3. The 5/i-Mixing Mechanism 

As described above, this mixing mechanism is activated when the hydrogen-burning 
shell approaches the homogenized region left behind by the retreating surface convection 
zone. 3 He is among the most fragile nuclei present, and the reaction 

3 He ( 3 He, 2p) 4 He , (1) 

has the unusual characteristic (among fusion reactions in stars) of lowering the mean molec- 
ular weight /x, creating a localised /x-inversion. This has already been noted by Ulrich (1972) 
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although he discussed it only in the context of burning during the pre-main-sequence con- 
traction phase. 

This /^-inversion is only about 0.025 R & above the hydrogen-exhausted core, and has 
the lowest molecular weight in the entire star (Fig 2). We emphasise that the 3 He-burning 
that we are discussing produces an inversion only because of the previous homogenisation 
of the outer layers by the FDU. In the absence of the FDU there would be no /z-inversion, 
because the 3 He peak would be superimposed on a 1 H/ 4 He gradient that gives a strongly 
stabilising /x-gradient. However when the entire outer layers are homogenised by the FDU, 
and then later encroached on by the gradual outward (in mass) motion of the burning shell, 
the background gradient has been removed and so the inversion can develop. 

Figure 2 shows the development of this /z-inversion as seen in a ID code, with no 
instructions for mixing in such regions. The maximum increase in A(l/yu) is seen as ~ 
5 x 10~ 4 . In the absence of 5/i-mixing, the magnitude of the inversion depends on the 3 He 
abundance. The temperature in the region where 3 He is burned is > 10 7 K, so that the 
material is fully ionised. The mean molecular weight is 

1 _ ^ (Z, + 1)X, 

The 3 He cross-section is larger than other rates, resulting in the production of 4 He and *H 
and reducing the molecular weight. Through reaction (1) fi decreases by an amount that 
depends on the change in the 3 He mass fraction. If the mass X 3 of 3 He in unit mass of gas 
increases by an amount 5X 3 (a negative quantity, in the present context) then according to 
equation (1) X\ increases by —2^X3/6 and X4 increases by —45X3/6, so that the change in 
H from equation (2) is given by 



1\ _ 2 25X3 35X3 3 45X3 5 ± = 5X3 

// ) ~ 16 1 (i /r ' 6 ' 

Figure 1 shows that the initial value of X 3 is about 2 x 10~ 3 and hence we predict a 
maximum A(l/yu) ~ 3 x 10~ 4 in good agreement with the 5 x 10~ 4 seen in Figure 2. It 
was this model, with the maximum possible value for A(l//x) that was mapped into the 3D 
Djehuty code. The resulting layer was dynamically unstable and began to rise at a relatively 
high speed ~ 10 2 m/s. 

A simple buoyancy argument was used to estimate the rise rate of these clouds to be 
of order 300 m/s, only modestly slower than the speeds expected in the convection region 
itself, and in agreement with the speeds seen in the 3D calculation. As these clouds rise 
(Figure 3), they are replaced with 3 He-rich material, and the process is continued. If such 
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speeds are maintained, the low-molecular-weight material will reach the convection zone in 
a few months. This is to be compared to the hundreds of millions of years required for a 
low-mass star to reach the helium flash after this ^//-mixing process begins. 

The disparity of the mixing to evolutionary timescale is such that Paper III used an 
instantaneous mixing approximation and the temperature/density structure of the stable 
region to estimate the amount of 3 He-processing. In the time taken by a 1 M & star to evolve 
up the giant branch, the 3 He in the envelope was expected to decline by about 3 e-folds or 
95%. A similar estimate for 12 C-processing indicated that 12 C would decrease by about 8%. 
While modest, this is sufficient to reduce the 12 C/ 13 C ratio from 25 to near 15. 

This analytic calculation suggested that the <5//-mixing might solve more than the conflict 
between the Big Bang nucleosynthesis and stellar evolution as to which produced the 3 He. It 
has the potential to solve many of the CNO isotope anomalies. Weiss & Charbonnel (2004) 
have identified the point where CNO composition differences (from expected values) begin 
to occur as 'where the hydrogen-burning shell encounters the deepest point to which the 
convective envelope ever reached', and this is just where the ^//-mixing mechanism becomes 
operational. 

In reality, however, the estimate in Paper II and the mixing seen in the 3D calculation 
is likely a substantial overestimate for the following reason. Once the 3 He starts to burn it 
will drive mixing of the low-// material and hence the actual //-inversion will be less than the 
maximum case discussed above. The actual value of 5fi will depend on the mixing speed, 
which we try to estimate below. 

To examine this phenomenon quantitatively, we have incorporated a ^//-mixing model 
into a ID stellar evolution code. The code tracks 16 isotopes, including all stable isotopes 
of the CNO tri-cycle: 1 H, 2 H, 3 He, 4 He, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 20 Ne, 21 Ne, 22 Ne, 
23 Na, x n. They are coupled through reaction rates taken from Caughlan & Fowler (1988). 
Following their recommendation, in reactions like 17 0(p,o;) 14 N and 17 0(p,7) 18 F that have 
an uncertain factor (0 to 1) on certain states, the factor was chosen to be 0.1. As discussed 
by Dearborn (1992), these factors are significant for the expected oxygen isotope ratios. In 
the sections below, this model will be tested and used to explore a range of masses and 
metallicities. Uncertainties in the reaction rates, and their effect on stellar evolution, have 
been discussed recently by Herwig et al (2006). 

We sum up the overall mixing that we expect as follows: 

(a) A stably stratified layer reaches temperatures where 3 He begins to fuse with itself, de- 
creasing the molecular weight of the layer in situ. Rather than develop as a dynamical 
instability (identified as a Rayleigh- Taylor instability in Paper II) the mixing is more like 
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a thermohaline process, with the buoyancy determined by the competition between the dif- 
fusion of heat from the layer on the one hand, and the mixing ("diffusion") of the 3 He fuel 
on the other hand. To the extent that the mixing is a diffusive process (and it is not) one 
could call this a "double diffusive" mixing, much like the thermohaline mixing seen in salty 
water (Stern 1960; Veronis 1965). In the case of salty water, the mixing is determined by 
the competition between the diffusion of salt on the one hand and heat on the other. In our 
case the competition is between the mixing of the material and the heat diffusion. 

(b) In addition, normal convection drives mixing throughout the classical surface convection 
zone. This mixing is much more rapid than the 5 fx- mixing discussed above. Classical con- 
vective mixing is likely to be reasonably well-modeled by a diffusion process, with a diffusion 
coefficient D ~ vol (Eggleton 1973, 1983) where vo is the mean turbulent speed of convection 
and I is a 'mixing- length' comparable to the pressure scale height. 

(c) Process (a), however, is clearly not a process that should be well-modeled by diffusion. It 
might be more reasonably considered advective rather than convective, and perhaps better 
characterised by a speed rather than a quantity like D of dimensions speed x length. This 
is because we expect finger-like structures to form, which will slowly mix and will not really 
resemble a homogeneous mixing process, as results from the diffusion equation. In the next 
two Sections we consider some estimates of the speed, and an estimate for an artificial D 
that will lead in practice to the sorts of speeds that we estimate. 

A number of papers have discussed mixing in the presence of a molecular weight inver- 
sion. By analogy with the situation seen in salt water, where surface evaporation creates a 
layer of hot salty water atop fresher cooler water, such situations are usually called "ther- 
mohaline" mixing. A theory for this mixing has been developed by Ulrich (1972) as well 
as Kippenhahn et al (1980) - see also the discussion in Kippenhahn (1974). The key result 
from these studies is that the mixing time-scale varies as the square of the size, d, of the 
unstable region. In the limit where d approaches the size of the star, the mixing timescale 
becomes the Kelvin-Helmholtz timescale. A rather unfortunate feature of the theoretical 
work is that estimates of the diffusion co-efficients vary by almost two orders of magnitude 
between authors, reflecting uncertainties in the geometry of the motion. Rather than place 
the emphasis on the existing attempts at modelling thermohaline mixing in stars, we prefer 
to take a phenomenological approach and estimate mixing speeds from first principles. 
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4. The Speed of 5/i-Mixing 

Our Djehuty calculation began from an artificial situation, in that the /z-inversion had 
been allowed to grow to its maximum size already in the previous ID calculation. We 
therefore saw rapid motion setting in rather quickly. In practice, the motion should have 
started as soon as the inversion began, and the motion set up by that, though slower, would 
have prevented the inversion building up to the size that we see in Figure 2. We attempt 
to make a somewhat more realistic estimate of the slower motion that we would expect to 
be set up in a roughly steady state. (Note that the 3D simulation did show rising blobs 
of material, but did not have the resolution to provide us with any guidance to the aspect 
ratio of the expected "fingers". We may investigate this in the future, although Djehuty is 
unlikely to provide much guidance on this question because of the expected long timescale 
being well beyond that of a hydrodynamic code like Djehuty. 

In a real stellar environment the buoyancy will begin when only a fraction of the 3 He 
is processed, and it will move only a short distance before coming to hydrostatic equilib- 
rium. When the bubble reaches this new equilibrium, its temperature will be lower than its 
surroundings by an amount 

^ = S ± (A) 

Heat will diffuse in, and the bubble will continue to rise on a thermal timescale. As a crude 
overestimate we can take the mixing time to be of the same order as the thermal timescale 
r, defined as the time it takes for the thermal energy in the radiative layer (3.65 x 10 46 ergs) 
to be replaced given the current luminosity (33.2L©), giving r ~ 5000 yrs. This gives an 
estimate (probably a considerable underestimate) of the mixing speed: 

Ar 

v > 0.2cm/s . (5) 

T 

where Ar is the thickness of the zone (~ 1 R Q ) between the /x-minimum and the conventional 
surface convection zone. While this mixing is not a diffusion process, the time for material to 
rise through the stable region is orders of magnitude shorter than the time for the hydrogen 
shell to burn through that material. As at result, this region maintains homogenization with 
the convective envelope. This allows a simple diffusion approach, and in the next section, we 
develop a diffusion approximation based on <5/i. We will show that increases in this diffusion 
coefficient by a factor of 100, or decreases by a factor of 30, have only a modest effect (~ 30%) 
on the determined abundances. 

We attempt to make a somewhat sharper estimate of the speed with which the bubble 
will rise. Consider a bubble of radius / containing lower-/i material. To obtain an upper 
limit to the velocity we will assume that the moving bubble is optically thin. The rate with 
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which heat energy enters the bubble can be estimated as 

F ~ 4tt/ 2 . acT 3 5T ~ 47r/ 2 acT 4 ^ . (6) 

As flux enters the low-temperature region, the temperature increases at a rate that depends 
on the volume: 

4 3 3pN kdT 

F ~ - 7T/ — , ' 

3 2 [i dt ' v ; 

with iVo being Avogadro's number and fc Boltzmann's constant, and so 

dT 2acT 4 

H^^kl M ■ (8) 
The temperature gradient that must be overcome is 

dT T 

= pV rad5 p. (9) 

This is accomplished as energy flows into the bubble resulting in a rise velocity of 

dT/dt 2acT 4 \Su\ 
v ~ — ■ — ~ - — - . (10) 

dT/dr glpV rad n 

The process of thermally-driven buoyancy described by equation (10) should leave the 
temperature gradient radiative, but we expect this to lead to elongated vertical structures 
resembling 'salt-fingers' (Wilson & Mayle 1988, Dalhed et al. 1999). As energy diffuses 
into the outer portion of the bubbles, the material rises, exposing the inner material. These 
slender structures have more surface area per unit volume than the bubble model used here, 
allowing energy to diffuse more rapidly into the low-// material, and this results in higher 
rise velocities than our estimate here. 

The factor is set at the deepest level of the mixing region where the heat from 

the H-burning shell drives the 3 He reaction. This gives 

M - ^ - lfidXs - -iV 0P 3Vi? nuc ^, (11) 



jj, 6 Qv dt 2 \ 3 / 6v 

where R nnc is the thermal average of reaction cross-section times speed. Then combining 
(15) and (16), 

2 acT 4 N R nuc fx ( X 3 



= . (12) 



Evaluating these equations in the region where 3 He is burned for a 1 M Q Pop I model, 
gives an initial speed in the bubble formation region of ~ 2 cm/s. This is an order of 
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magnitude above the first crude (under-)estimate, but well within the range of diffusion 
coefficients tested in the next section. The speed in the warm processing region is limited 
by the conversion rate of 3 He and it is here that S/j, is established, following which it can be 
held constant and equation (10) used to calculate the speed. This speed will be typically 
somewhat larger than the initial speed. 

Averaged over the classically stable region, equation (10) in Pop I models show speeds 
of ~ 1 — 2m/s. These high speeds are irrelevant to the processing region and only serve to 
speed homogenization with the convection zone. Pop II models develop lower \8 [/,{/[/, values, 
and average velocities nearer 0.5 m/s. This estimate was repeated for several core masses 
as the model evolved up the FGB, resulting in mixing times of between 10 and 50 years. 
Models of different mass were also examined, with the result that lower-mass models mix 
somewhat more slowly and higher-mass models somewhat faster. Still, for Pop I models the 
mixing timescale was under 100 years (Figure 4). We note that the difference in mixing speed 
at the bottom of the envelope, and the exact location of the peak in the \x inversion, will 
determine the final carbon isotope ratio. As we discuss below, these details require further 
investigation and are beyond the scope of the current paper. 

Extreme Pop II models (Z = 0.0001) were also examined. The initial mixing in these 
models starts at higher temperatures, resulting in greater values of \5 [/,]/[/,. The mixing 
timescales start near 10 years for all the masses examined. However, as these models evolved 
the timescale rapidly increased to ~ 100 years. This is probably a result of the rapid de- 
struction of 3 He. A mixing time of 100 years is about 0.01% of the time to burn through the 
classically stable region. We can argue that the classically stable region will be homogenised 
with the classically convective surface zone on a short timescale. 

We do not feel that simplistic analyses such as equations (6) to (12) are a very reliable 
guide to the velocities to be expected, and in particular we do not feel that a roughly 
mixing-length-like theory is appropriate for the rate of travel of salt-finger-like mixing. In 
the process described after equation (10), based on 3D simulations in the papers quoted, the 
bubble does not in fact wait for all the heat to diffuse in to the center before rising starts. 
That appears to be a key factor in the finger-like character of the structures. Because of that 
we have used an optically thin approximation in equation (6) rather than an optically thick 
approximation. But several other steps in this kind of back-of-the-envelope estimation are 
quite problematic. We reiterate that although theories for thermohaline mixing exist, they 
are uncertain by typically two orders of magnitude, as reflected in the preferred diffusion 
co-efficients used by Ulrich (1972) and Kippenhahn et al (1980). 

We prefer (next Section) to use an alternative simplistic procedure, which is to choose 
a 5/i-mixing diffusion coefficient such that it will give speeds of advance comparable to the 
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kinds of speeds we have estimated above. We are fortunate that it does not matter, to within 
several orders of magnitude, just how rapid is the actual motion. As long as the motion is 
rapid compared to the very slow overall nuclear timescale of the star then the results will 
be essentially the same, and independent of the mixing speed. Whether the speed of rising 
elements is 2 m/s or 0.2 cm/s, this is still faster than the speed with which material would 
have to travel to mix on a nuclear timescale, ~ 10~ 3 cm/s. Indeed, this is confirmed below 
where we show that varying the mixing speed by more than three orders of magnitude 
provides quantitative changes of order ~ 30% in the determined surface abundances. 



5. The 5/i-Mixing Model 

In ID codes, mixing always requires some physical model to approximate the process. 
Our ID code treats normal convection as a diffusion process and solves a second-order 
equation for each isotope: 

dX\ d 2 , 2 dX n dX fdm\ , , 

■w)r^ D{4 * rp) ^ + R+ a^{w) t - (13) 

The first term in the equation for the rate of change of the isotope is for the convective 
diffusion, the next term (R) incorporates the nuclear reaction rates, and there is a final 
term to deal with the mesh motion, because the mesh is non-Lagrangian; k is the meshpoint 
number. 

In a standard convective region the diffusion coefficient D that we use takes a form 
that scales quadratically with the temperature-gradient excess over the adiabatic value, and 
inversely with the nuclear timescale: 

£ = f^[max(0,V r -V a )f , (14) 

''nuclear 

where V r and V a are the usual radiative and adiabatic temperature gradients from the 
mixing-length theory of convection. The timescale tnuciear is an estimate of the nuclear 
evolution timescale: t nuc i car = O.IEXqM/L, where E is the nuclear energy available from 
hydrogen burning, X is the abundance of hydrogen in the outer layers (~ 0.7) and M is the 
mass of the star: t nuc i C ar ~ 10 10 yrs for the ZAMS sun, and is ~ 10 s yrs in the middle of 
the FGB. The dimensionless and largely arbitrary factor -F conv is simply chosen to be a large 
number such that the composition in a convective region homogenizes in a time much shorter 
than the nuclear time scale (t nuc icar)- We find that a value of 10 6 works well. A larger value 
would be more physical, but would lead to numerical difficulty since the stepwise difference in 
composition becomes so small as to be susceptible to rounding error. There is little difference 
in practice between a convective region being homogenised in 10~ 6 t nuc i ea r and 10~ 10 t nuc i ear . 
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To model our 5/i-mixing process we have created a diffusion coefficient wherever there 
is an inversion. The form that we have used is simple, and in analogy with the previous 
equation: 

F inv r 2 

D = - ^ min ) (k > k mhl ) (15) 

^nuclear 

= (k < k min ) , (16) 

where (a) /i m i n is the smallest value of jj, in the current model, (b) /i m i n occurs at meshpoint 
number fc min , and (c) k, the meshpoint number, is counted outwards from the center. Thus 
we introduce mixing wherever there is an inversion in /z, in the same form as usually used for 
mixing, but with a factor Fi nv to be determined. In a region where both of equations (14) and 
(15) give D > 0, we only use (14). Again, the factor -Fi nv , if large, simply assures homogeneity 
in a time much shorter than the nuclear time scale. It remains for us to determine a value 
of F inv that is appropriate for 5/x-mixing. 

To estimate the speed that corresponds to a chosen value of F inv , we performed a 
numerical test in which a step function was installed in an element that was not being used 
in the nucleosynthesis network. The position of the step was located just above the point 
where the /z-inversion would form. Below this point, the mass fraction was set to 10~ 7 , and 
above this point dropped to 10~ 10 . Once the inversion develops, mixing begins, and a stable 
gradient is formed. The rate at which the material below the step is transported to the 
surface was then monitored to obtain an effective speed (with the usual convective diffusion 
coefficient turned off to avoid confusion). This is illustrated in Figure 5, for F inv = 10 2 . Some 
material is seen to have travelled ~ 6 x 10 11 cm in 80,000 yrs, ie. with a speed of ~ 0.3cm/s. 
This is about as slow as our slowest estimate above, and hence is chosen as our standard 
value. Increasing Fi nv by 10 times results in a speed that approached 1.5 cm/s. Increasing 
this arbitrary value by another factor of 10 (100 times the standard value) increases the speed 
to 6 cm/s. These speeds are sufficient to mix the stable region with the outer convective 
region on a time that is short in comparison to the evolution. These speeds are near the 
startup velocity found where the reduced | <5/x | / is established, where the nucleosynthesis is 
done, and where the composition changes are most sensitive to mixing speed. 

Figure 6 shows the diffusion coefficient that we used for <5/i-mixing, along with the 
diffusion coefficient used for the ordinary convective envelope. Beyond the point where the 
switchover occurs, we plot (and use) only our value for normal convective mixing. We note 
that 5/i-mixing produces diffusion co-efficients that are up to two orders of magnitude larger 
than found for rotating models (Palacios et al 2006). 

To illustrate the effect of this 5/i-mixing on a star's evolution, a 1 M (Z = 0.02) model 
was evolved from a pre-main-sequence configuration to the helium core flash. The initial 
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12 C/ 13 C ratio was chosen to be 90, and as the model reached the giant branch (M core ~ 
0.2 Mq) the first dredge-up reduced the ratio to 29.5. When the new mixing mechanism 
began, the 12 C/ 13 C ratio declined rapidly to near 15 (Figure 7). At this point, the reduced 
3 He abundance slowed the mixing to the extent that between a core mass of 0.3 and 0.45 M® 
the ratio declined only to 14.3. Nearly 93% of the initial 3 He was destroyed by the tip of the 
Giant Branch. 

When the 5/z-mixing is active, the striking molecular weight inversion (red in Figure 8) 
does not develop. Instead, a much more gradual profile is developed with a modest cusp 
in the region where the 3 He burning occurs. We find A(/z//i) ~ 10~ 5 but the details will 
depend on the assumed mixing speed as well as the temperature of the approaching shell. 



6. Sensitivity of the 12 C/ 13 C ratio and 3 He destruction to F inv 

It is important to realise that our preferred value of F inv was chosen to match our 
estimates of the mixing velocity. Nevertheless, we now wish to test the sensitivity of the 3 He 
destruction and the final 12 C/ 13 C ratio to the value of -F; nv . A series of runs were made on a 
1 M e Z = 0.02 model. The value of F inv was varied over a factor of 10,000 (from 0.01 of the 
standard value to 100 times the standard value). To test the sensitivity to mesh resolution, 
this test was done for models with both 300 zones and 750 zones. 

Table 1 shows the effect of varying F inv by factors of up to 10 4 . For values of F inv < 10 
the mixing is too slow and the 3 He burning is incomplete. For rates near the standard value 
chosen, destruction is near 90%, and 12 C/ 13 C ratios are a minimum. Larger diffusion coef- 
ficients only modestly reduce the 3 He abundance and give a little less reduction in 12 C/ 13 C. 
The differences are not significant. Note that the preferred value of F inv = 100 was chosen 
becasue it gives a mixing speed close to the minimum that we estimated; hence F[ nv = 10 
is an order of magntude slower than our estimate of the minimum. The table also shows 
a modest sensitivity to mesh resolution. This information is repeated more graphically in 
Figure 9. 

Next we examined a range of masses (for Z = 0.02, and 300 zones) comparing diffusion 
coefficients that differ by a factor of 100. For the higher diffusion coefficients the 3 He de- 
struction is almost identical except for the higher mass where less 3 He is produced in the first 
place. There are modest differences in 12 C/ 13 C ratios that at present should be considered 
as uncertainty in the modeling of the mechanism. 

Table 2 shows that the behavior seen in the 1 M Q model holds over the range of inter- 
esting masses. Between 0.8 and 2.0 M , higher diffusion coefficients result in small changes 
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in the 3 He destruction, and very modest differences in 12 C/ 13 C ratios. The differences in 
12 C/ 13 C ratios caused by varying mesh and mixing coefficient (±2 for Pop I models) should 
be considered uncertainty in the modeling of the mechanism (at present). When the ratio 
drops to near the equilibrium value of 3.5, as we will see in Pop II models, these factors have 
much less effect on the 12 C/ 13 C ratios. 



7. 12 C/ 13 C and 5/i-Mixing 

The following section examines the 12 C/ 13 C ratios, and the helium production, for a 
range of low-mass stars with Pop I and Pop II metallicities. We use the standard diffusion 
coefficient (F inv = 100) developed above for the <5yU-mixing. Consistent with Anders & 
Grevesse (1989), the initial 3 He mass fraction was taken to be 2x 10~ 5 . While it is appropriate 
to use a higher value to account for the conversion of D to 3 He on the pre-main sequence, 
this difference is minor when compared to the main sequence production of low-mass stars. 

Table 3 shows the 12 C/ 13 C ratios for models ranging from 0.8 to 2.0 M , and for metal- 
licities from solar to l/50th solar. In the absence of an additional mixing process, the final 
(tip of Giant Branch) value of the carbon isotope ratio depends on mass. As the mass rises 
from 0.8 to 2.0 M Q , the expected 12 C/ 13 C ratio drops from near 35 to near 20, with very little 
dependence on Z. This mass dependence is seen in the FDU columns of Table 3, showing 
the 12 C/ 13 C values after FDU and before 5/i-mixing begins. Once the mixing begins, the 
12 C/ 13 C ratio rapidly drops to a lower value, and the final range of ratios ('Final' in Table 3) 
show a considerably reduced range. For solar metallicities, stars in this mass range all show 
12 C/ 13 C « 14.5 ± 1.5. Similarly, for Z = l/50th solar, 12 C/ 13 C « 4.0 ± 0.5. The 2.0 M Q 
models are not included in these averages, as the 5/i-mixing begins just prior to helium core 
flash in the Pop I model, and has not begun in the Z = 0.0004 model. 

That the final 12 C/ 13 C values converge for a broad range of masses is an interesting 
result, and is shown graphically in Figure 10. Before 5/i-mixing begins, the carbon isotope 
ratios show the usual mass dependent range, but not afterwards. To illustrate Z dependence 
in more detail, a star of mass 0.9 M & was evolved with various values of Z between solar 
and l/200th solar (Figure 11). The 12 C/ 13 C ratio is seen to vary smoothly from 14.8 to 3.5. 
Additionally there are big changes in 14 N/ 15 N, and small changes in O isotope ratios (mostly 
due to 18 and 17 0). Table 4 provides the same information in tabular form. 

In the low-mass stars of interest here, the PP chain dominates evolution on the main 
sequence, but the hydrogen-burning shell on the giant branch operates on the CNO cycle. 
Stars with a lower metallicity Z will also have fewer CNO nuclei, so that the shell must 
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burn at a somewhat higher temperature for the same energy production rate. Additionally, 
the penetration of the surface convection region is not as deep at low Z, and the inversion 
is not initiated until the core grows to a larger mass. These effects combine to result in 
a higher temperature in the place where the molecular weight inversion develops. The 
position of this inversion is the result of competition between the advance of the hydrogen- 
burning shell (whose temperature is dependent on the star's metalhcity) and the speed of 
mixing. As Z decreases from 0.02 to 0.0001, the core mass at which the inversion occurs 
increases (from 0.233 to 0.358 for models of mass 0.9 M ). The luminosity, logL, at the 
start of mixing increases from 1.4 to 2.4, and the temperature at the base of the mixing 
increases from 16.5 million K to 24.0 million. As a result, the mixing begins at a metallicity- 
dependent temperature, and the 12 C/ 13 C ratio achieves different values before the reduced 
3 He abundance slows the process (Figure 12). 

It is again important to note that the composition dependence of our results was in no 
way specified by us, but is a direct result of the composition dependence of the temperature 
in the burning regions. All we have done is determine a mixing velocity (with no explicit 
composition dependence). The result is exactly as observed: the carbon isotope ratios fall 
further for lower metallicities. 

As a final note, these models were all run with a mixing length that was 1.8 pressure 
scale heights (to fit the solar radius). Changing the mixing length ± 0.2 resulted in changes 
in the 12 C/ 13 C ratio that were <0.4. 

8. The Metallicity Effect and Observations 

Gilroy & Brown (1991) measured the carbon isotope ratios in the stars of M67 (near 
1.2 Mq). They report that 'the subgiants seem to have undergone little or no mixing', and 
that the lower giant branch stars exhibit 'normal first dredge-up mixing ratios'. However 
they find 8 upper giant branch and clump giants with 12 C/ 13 C ratios between 11 and 15. This 
is in excellent agreement with the value near 13.5 that results from 5/i-mixing (Table 3). At 
the lower metal enrichments, Pavlenko et al. (2003) observed giants in the globular clusters 
M3, M5, and M13 ([M/H] = -1.3, -1.4, and -1.6 or Z < 0.001), finding 12 C/ 13 C between 3 
and 5. Again this is in excellent accord with our models (Table 3). They also observed M71 
([M/H] = —0.71 or Z « 0.004) finding less processing. Here the 12 C/ 13 C ratios show values 
near 7 (5 to 9). The 0.9 M model that we used to illustrate the metallicity effect gives an 
expectation for 12 C/ 13 C ratios near 8. 
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9. 3 He and ^-Mixing 

As stated earlier, in the calculations done here the initial mass fraction of 3 He was 
taken to be 2 x 10~ 5 . In Table 5, a 0.9 M & model is evaluated for a range of metallicities 
to find the peak enhancement in the mass of 3 He, M( 3 He, peak)/M( 3 He, init). For solar 
composition, the enhancement reaches a factor of 61. With shorter main sequence lifetimes, 
Pop II models of this mass reach only about 38. Increasing the main sequence value of 3 He 
by an order of magnitude, to 2 x 1CT 4 , by assuming all of the deuterium is burned to 3 He on 
the pre-main sequence makes surprisingly little difference. For Pop I abundances, starting 
with an enhanced initial 3 He abundance associated with converting 2 H to 3 He results in a 
peak enhancement of 64 times the ISM 3 He used to form the star (2 x 10~ 5 ) instead of 61. 
For an extreme Pop II metallicity, the peak enhancement is 41 instead of 38. 

These very large peak enhancements were problematic in reconciling stellar nucleosyn- 
thesis and Big Bang nucleosynthesis with observed abundances. With <5/i-mixing the final 
enhancement for Pop I abundances is reduced to 3.1 (3.4 when 2 H conversion is included), 
and for Pop II models the final value is 0.52 (0.53 with 2 H conversion). 

Table 6 shows the same behavior for stars in the mass range of 0.8 to 2M . Table 6 
shows the ratio of the peak Pre-mix mass of 3 He to the original mass of 3 He in the star, as 
well as the post-mixing ratio 

As found in earlier papers, the greatest potential 3 He enhancement occurs for stars of 
1 M & and below. In these models, it is the large 3 He enhancement that enables the <5/z-mixing 
to destroy 90 to 95% of the potential 3 He contribution to the interstellar medium. As a final 
note, at and above ~ 2 M & the 5/i-mixing operates to (near) completion only for the Pop I 
model. In the Z = 0.001 calculation, the mechanism begins just prior to the helium flash 
and is incomplete. In the Z = 0.0004 model, the helium core flash precedes any significant 
^//-mixing. 



10. 10 O and 23 Na with <5/i-Mixing 

Some observational evidence has suggested that 16 and 23 Na may change on the upper 
giant branch. In the 1.5 M models run with Z = 0.02 and 0.0001, there was substantial 
23 Na enhancement: <5( 23 Na) 19% and 68%. However, this was just the usual enhancement 
expected from main sequence processing followed by giant branch mixing, and is not signifi- 
cantly effected by 5/x-mixing. When we examined the 0.9 M Q model the greatest change was 
seen for the extreme Pop II abundance (l/200th solar), where the 23 Na abundance climbs 
by 3%. Only about half of this change is due to 5/i-mixing. The 16 depletion is a trivial 
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0.25% (Figure 13). 

In low-mass, low-metallicity stars, the 5/i-mixing brings the effective bottom of the 
convection zone to a point where 16 and 23 Na are beginning to change. Because of this the 
change in these isotopes is dependent on the location that the mixing model finds for the base 
of the 5/i-mixing, where the mixing and burning of 3 He are in balance. The mixing model 
developed here does not appear sufficient to explain an observable 16 depletion or 23 Na 
enhancement, but this result warrants additional investigation. Because of the temperature 
sensitivity, a modest amount of (downwards) overshoot, or turbulent mixing, could change 
this result. Indeed, we note that our 5/i-mixing mechanism makes it far easier for other 
mixing mechanisms to have an effect. In the absence of 5/i-mixing there is a radiative region 
of approximately 1 R Q (or 0.02 M or 10 pressure scale-heights) between the hydrogen shell 
and the convective envelope. Any proposed mixing mechanism must lift material form the 
shell through this radiative region to make contact with the convective envelope. However, 
with 5/i-mixing acting from just above the shell, and linking to the convective envelope, the 
radiative region between the shell and the bottom of the mixed region is now reduced to 
approximately 0.005 R Q (or 0.0004 M Q or less than one pressure scale height). Thus other 
forms of mixing may be aided by the operation of 5/i-mixing. 



11. Mass Loss after the Mixing does not Matter 

As noted in the introduction, the blue horizontal branches of globular clusters indicate 
that a substantial fraction of the envelope is lost prior to the helium core flash. This was 
the origin of the excess 3 He production by these stars. With this new mixing mechanism, 
once the /z-inversion occurs, the 3 He abundance drops rapidly. As the mixing depends on the 
square of the 3 He abundance, the process slows dramatically when 90% of this isotope has 
been consumed. Mass loss after this rapid drop does not cause contamination problems with 
excess 3 He, and leads to little change in the surface abundances or yields. To illustrate this, 
a 0.8 M Q Z = 0.0001 model was run with no mass loss, and then again with 0.2 M of mass 
loss near the tip off the giant branch (and after the <5/i-mixing) . The surface abundances 
show little change from the mass loss. 

In the absence of <5/z-mixing, the mass losing model ejects 52 times as much 3 He into 
the ISM as it took when it formed. This was the basis of the problem of reconciling Big 
Bang nucleosynthesis with yields from stellar evolution. With 5/i-mixing, the mass loss ejects 
l/3rd of the 3 He into the ISM as it took when the star formed, and retains only l/7th of the 
intake to be further processed or ejected. 
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12. Conclusions 

Our first conclusion is that <5/z-mixing is a significant and inevitable process in low- 
mass stars ascending the giant branch for the first time. Once it begins, the timescale is 
short (compared with the overall nuclear timescale of ~ 10 8 yrs, and it maintains a nearly 
homogeneous composition down to base temperatures in the region of 16 to 25 million K, 
allowing nuclear processing. The result is an observable change in the expected abundances 
of 3 He and the CNO isotopes. 

This mixing mechanism is driven by the destruction of 3 He, and is self-limiting. The 
lowest mass stars (< 1.25 M ), that were expected to produce a problematic excess of 3 He, 
quickly destroy 90 to 95% of that isotope. As a result, the 3 He returned to the ISM is within 
the limits posed by Hata et al (1995). This mixing also modifies the 12 C/ 13 C ratios. Instead 
of showing a significant mass dependence in that ratio, we find a metallicity dependence 
instead. Shortly after this mixing begins, Pop I stars between 0.8 and 2.0 M should all 
drop to a value near 14.5. Extreme Pop II stars in this mass range should show ratios near 
4. Likewise the nitrogen and oxygen isotope ratios are substantially affected (Figure 14). 

As a first effort to develop a ID model for this mixing process, we have tried to validate 
it with reference to the original 3D modeling, and basic physical arguments. At the moment, 
variations due to mesh resolution and mixing speed suggest an uncertainty of ±2 for Pop I 
values of the 12 C/ 13 C ratio (much less for Pop II, where it reaches its equilibrium value of 
~ 4). The fractional destruction of 3 He seems less sensitive to such choices. We have also 
held the initial isotope ratios fixed ( 12 C/ 13 C = 90, 14 N/ 15 N = 270, 16 0/ 17 = 2625, and 
16 0/ 18 = 490), and different values must be run to study chemical evolution. 

Although our diffusive model of <5/i-mixing is very tentative, we have endeavored to 
show that the details are not very important. We estimate mixing speeds between 0.2 and 
2 cm/s in the burning region, increasing to perhaps 1 m/sec in the radiative zone below the 
convective envelope. We also showed that any rate within one or two orders of magnitude of 
the rate we used is going to produce much the same result so far as mixing is concerned. It 
is rather unsatisfying not to have a more accurate estimate of the speed of mixing. However, 
for the problems discussed in this paper the speed is not crucial, as long as it is fast relative 
to the nuclear timescale, and certainly this is what we find, by some orders of magnitude. 
For other applications of the same mechanism it may be that a more accurate estimate of 
the speed is required. We note that appealing to current theories of thermohaline mixing 
is unlikely to help, as these vary from one author to the other by two orders of magnitude. 
Further, since the mixing timescale is essentially the Kelvin-Helmholtz timescale, then a 
direct numerical simulation is impossible at present. 
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We believe that our model accounts remarkably well for at least two results: the fact 
that 3 He is not as much enhanced in the ISM as earlier models suggested, and the fact that 
the 12 C/ 13 C is seen to decrease substantially further than was found in earlier models. We 
expect in a later paper to look in more detail at oxygen isotopes as well as 7 Li in red-giants. 

Finally, this process does not appear to have significant impact on the 16 depletion 
or 23 Na enhancement, but with the temperature sensitivity (and rate uncertainties) of these 
rates, it comes close. Any overshoot in the mixing at the bottom of the region will be 
important. Also, the thickness of the stable region that surrounds the hydrogen burning 
shell is reduced from over a solar radius to a few hundredths of a solar radius. This may 
enable other mechanisms (rotation, magnetic fields, ...) to create variation in the observed 
abundances. Alternatively, the homogeneity seen in clusters like M67 (Gilroy & Brown 1991) 
might be used to limit models for such mixing. 
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Fig. 1. — Composition structure of a 0.85 M Q Pop I model with X = 0.7, Z = 0.02. Left 
panel - near the end of its main-sequence life. Right panel - near the deepest penetration 
of the surface convective region on the giant branch. X H - vermilion; 4 He - magenta; 12 C - 
black; 16 - blue; 14 N - green; 3 He - red dots; 13 C - black dots. 
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Fig. 2. — The development of a molecular-weight inversion is shown (versus mass) for a 
1.0 M® Pop I model with X = 0.7, Z = 0.02. curves are shown for million-year 

increments beginning at 11.67 Gyr. 



Fig. 3. — The location of the hydrogen burning shell is shown by a blue contour at a fixed 
mass fraction of 13 N. Shown in red are rising clouds in which the hydrogen abundance is 
marginally higher than in the surrounding material. 
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Fig. 4. — Mixing times along the giant branch, for Pop I and Pop II models with a range of 
masses. 



-26- 




Fig. 5. — The mass fraction profile created by the proposed diffusion coefficient in a 1 M & 
star [Z = 0.02), starting from a pure step-function profile. Each curve is separated by 10,000 
years. For the value of F inv selected here (100), the speed at which the material is carried 
outward corresponds to about 0.3 cm/s. 
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Fig. 6. — The diffusion co-efficient that we use (in cm 2 /s), (a) solid line, cycle 741, before our 
mechanism sets in, so that there is only the normal surface convection zone (D determined 
by equation 14); (b) dotted, cycle 752, just after 5/i-mixing begins, D from equation (15) 
except for the convective region, which uses equation (14); (c) dashed, cycle 821, well after 
the mixing is established. Radius is in solar units. 
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Fig. 7. — The 12 C/ 13 C ratio versus core mass in a 1 M model. 
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Fig. 8. — Approximately 3 million years after the hydrogen burning shell approaches the 
homogeneous region, the profiles for a 1M model with no mixing (red) and a 5fi- 
mixing model (blue). 



-30- 



30- 



25 



20- 



u 



10- 



5 - 



0^ 



_i i i i i i i i i i i i i i i i i 



Speed range of 
0.2 to 2.0 cm/s 

300 Zones 

750 Zones 




0.0001 0.001 0.01 0.1 1 

Multiplier on F jnv 



10 



100 



Fig. 9. — Surface carbon isotope ratios as a function of F inv and the number of zones used 
in the ID model. The region corresponding to the estimated speed range of 0.2 to 2 cm/sec 
is shown by the vertical green lines. The x-axis is the multiplier used on the standard 
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Fig. 10. — 12 C/ 13 C ratio for various masses using our standard F inv = 100. The nearly 
vertical dark line marks the onset of the 5/x-mixing process. 
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Fig. 11. — Post-mixing, surface isotope ratios in a 0.9 M Q model with a large range of initial 
abundances. 



- 33 - 




-34- 




Fig. 13. — Change in the surface abundance of 16 and 23 Na in a 0.9 M model evolved with 
various metallicities. 
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Fig. 14. — CNO isotope ratios in a 1M (Z = 0.0001) model are plotted against core mass. 
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Table 1. Effect of Varying F inv and Number of Zones 



300 Zones 750 Zones 

F inv 12 C/ 13 C 3 He Destroyed 12 C/ 13 C 3 He Destroyed 






29.5 


0% 


29.1 


0% 


1 


18.9 


59.0% 


27.8 


35.8% 


3 


14.0 


75.8% 


22.9 


46.7% 


10 


12.2 


87.3% 


16.1 


64.9% 


30 


12.3 


91.3% 


12.2 


80.6% 


100 


13.4 


92.6% 


11.1 


89.1% 


300 


15.7 


92.4% 


11.7 


91.9% 


10 3 


16.1 


92.8% 


12.7 


93.1% 


10 4 


16.5 


93.9% 


14.1 


93.8% 
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Table 2. Effect of varying mass and Fj 



Mass FDU 12 C/ 13 C 3 He Destroyed 12 C/ 13 C 3 He Destroyed 



0.80 36.9 15.9 96.4% 19.5 96.7% 

0.85 34.0 15.3 95.7% 18.5 96.0% 

0.90 32.2 14.5 94.8% 17.6 95.4% 

1.00 29.5 13.4 92.6% 16.5 93.5% 

1.25 25.6 13.0 85.7% 14.9 89.1% 

1.50 23.6 13.7 74.9% 14.4 82.3% 

2.00 22.3 17.0 45.1% 14.9 63.6% 
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Table 3. 12 C/ 13 C ratios 





X=0.70, 


Z=0.02 


X=0.738, 


Z=0.001 


X=0.74, 


Z=0.0004 


Mass 


FDU 


Final 


FDU 


Final 


FDU 


Final 


0.80 


36.9 


15.9 


34.1 


5.3 


35.0 


4.2 


0.85 


34.0 


15.3 


31.5 


5.0 


31.8 


4.0 


0.90 


32.2 


14.5 


29.6 


4.9 


30.0 


4.0 


1.00 


29.5 


13.4 


27.3 


4.9 


27.4 


4.0 


1.25 


25.6 


13.0 


24.3 


5.0 


24.3 


4.1 


1.50 


23.6 


13.7 


24.3 


5.2 


22.7 


4.6 


2.00 


22.3 


17.0 


21.2 


14.2 


21.0 


21.0 
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Table 4. Surface Isotope Ratios at the Top of the Giant Branch 



z 


3 He/ 4 He 


12 C/ 13 C 


14 N /15 N 


16Q/17Q 


16Q/18Q 


0.02 


1.8xl0 3 


14.5 


0.8xl0 4 


2.5xl0 3 


513. 


0.01 


2.3xl0 3 


11.3 


1.4xl0 4 


2.4xl0 3 


521. 


0.004 


3.2xl0 3 


8.0 


1.8xl0 4 


2.3xl0 3 


537. 


0.001 


5.2xl0 3 


4.9 


2.1xl0 4 


1.9xl0 3 


584. 


0.0004 


6.8xl0 3 


4.0 


2.2xl0 4 


1.7xl0 3 


640. 


0.0001 


9.1xl0 3 


3.5 


2.3xl0 4 


1.5xl0 3 


752. 



-40- 



Table 5. 3 He production in a 0.9 M model. 



z 


0.02 


0.01 


0.004 


0.001 


0.0004 


0.0001 


Peak 


61.1 


56.5 


48.7 


40.0 


38.7 


38.1 


Final 


3.1 


2.3 


1.6 


1.0 


0.7 


0.5 


Change 


94.8% 


95.8% 


96.7% 


97.6% 


98.1% 


98.6% 
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Table 6. 3 He/ 3 He (original) 





X=0.70, 


Z=0.02 


X=0.738, 


Z=0.001 


X=0.74, 


Z=0.0004 


M 


Peak 


Mixed 


Peak 


Mixed 


Peak 


Mixed 


0.80 


76.6 


2.7 


54.7 


0.86 


53.0 


0.62 


0.85 


68.2 


2.9 


46.5 


0.90 


45.1 


0.66 


0.90 


61.1 


3.1 


40.0 


0.95 


38.7 


0.73 


1.00 


49.7 


3.7 


32.9 


1.12 


30.0 


0.83 


1.25 


31.7 


4.5 


21.8 


1.60 


19.9 


1.21 


1.50 


21.8 


5.5 


15.7 


2.08 


14.4 


1.91 


2.00 


12.8 


7.0 


9.5 


6.70 


8.7 


8.39 
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Table 7. Enhancements with and without Mass Loss 



Isotope Ratios 


No Loss 


Loss 


3 He/ 4 He 


9065.5 


9466.7 


12 C/ 13 C 


3.5 


3.5 


14 N/ 15 N 


22564.1 


22330. 


16Q/17Q 


1913.4 


1920.8 


16Q/18Q 


702.6 


677.8 


5( 16 0) 


-0.01% 


-0.01% 


5( 23 Na) 


1.2% 


1.3% 



